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Abstract 



Cosmic superstring loops generically contain strings of different tensions that meet at 
Y-junctions. These loops evolve non-periodically in time, and have cusps and kinks that 
(N) ■ interact with the junctions. We study the effect of junctions on the gravitational wave 

signal emanating from cosmic string cusps and kinks. We find that earlier results on the 
strength of individual bursts from cusps and kinks on strings without junctions remain 
largely unchanged, but junctions give rise to additional contributions to the gravitational 
wave signal coming from strings expanding at the speed of light at a junction and kinks 
r~| \ passing through a junction. 

<-H ■ . 

> • 1 Introduction 

! In string theory models of brane inflation, fundamental (F) and D-strings are produced at 
the end of inflation when the branes collide and annihilate [U [21 El H]. In certain scenarios 
the resulting superstring networks are stable and expand to cosmic size, with predicted string 
tensions in the range 1CT 11 < Gfi < 10~ 6 . This raises the possibility that cosmic superstrings 
could provide an observational signature of string theory. Indeed it has been argued that the 
gravitational wave (GW) signals from cusps of oscillating loops on cosmic strings should be 
^ | detectable by LIGO and LISA for string tensions as low as 10~ 13 [1J EJ Ej- 

However the GW predictions of [B] may not be directly applicable to cosmic superstring 
networks, because F and D-strings form an interconnected network in which they join and 
separate at Y-junctions. Each junction joins an F-string, a D-string and their bound state. 
This means that closed loops containing junctions in a cosmic superstring network do not 
evolve periodically in time. Furthermore kinks interact with junctions, which leads to novel 
contributions to the GW signal as well as a more complicated network evolution. For these 
reasons we re-examine here the gravitational wave signal emanating from cusps and kinks on 
cosmic string loops taking in account the presence of junctions. We note that in a similar 
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spirit, Damour and Vilenkin have calculated the effect on GW signatures of cosmic strings of 
a reduced reconnection probability of intersecting strings and of a reduced typical length of 
newly formed loops. It was found that earlier results obtained for field theory cosmic strings 
remain largely valid for a rather wide range of network parameters. 

We first set up the formalism for calculating GW bursts emitted by cusps and kinks on 
non-periodic cosmic string loops. The non-periodic evolution of loops with junctions renders 
the GW calculation somewhat more involved because one can no longer factorize the Fourier 
transform of the GW amplitude. We show one can nevertheless integrate the stress energy 
over the string worldsheet and obtain an analytic expression for the high frequency behavior 
of the various contributions to the GW bursts. We find that earlier results on the strength of 
individual bursts from cusps and kinks on strings without junctions remain largely unchanged, 
but junctions give rise to additional contributions coming from strings expanding at the speed 
of light at a junction and kinks passing through a junction. We analyze the latter contributions, 
which provide a possible observational discriminant between ordinary cosmic strings and cosmic 
superstrings. 

We concentrate here on the calculation of individual gravitational bursts. This has the 
advantage of being a fully tractable exercise. By contrast, the observable signal of such bursts at 
the present time depends also on the details of the loop evolution as well as on the cosmological 
evolution of the network. In particular junctions tend to enhance the number of kinks on loops. 
One might expect this amplifies the GW signal from cosmic superstring networks compared 
to the signal from cosmic string networks without junctions. A more detailed analysis of this 
effect will appear elsewhere [T2] . 

2 GW emission from cosmic superstrings with junctions 

In this section we consider the GW emission from a cosmic superstring loop, and generalize the 
results of Damour and Vilenkin (5j El [7] on the emission of gravitational bursts from standard 
loops without junctions. 

We consider closed loops of cosmic superstrings containing two Y-j unctions, as shown in 
figure [TJ (Generalization to loops with 4 or more junctions is straightforward in principle). The 
loop consists of three local Nambu-Goto cosmic strings of tensions \i q (q = 1, 2, 3), which meet 
at two Y-junctions where the junctions are labeled by A and B. The action describing the 




Figure 1: Loop formed by 3 strings and 2 junctions 
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system has been set up and analysed in [8], and is given by 



S = ~ E IH [ dt ^ d ^ \J< (1 - ^) + E E / dtf/-[x,(t,s g J (t))-X J (t)]. (1) 

5=1,2,3 7 J=(A,S) 9 7 

The first term is the Nambu-Goto action for each of the three strings in the loop, where we have 
assumed a flat spacetime geometry with signature ( — h ++) and used the standard conformal- 
temporal gauge to parametrize each string's worldsheet. Thus each string is described by its 
spatial coordinates x q (a q ,t) where t coincides with background time. (As discussed in jT3], 
cosmic superstrings should be described by the Dirac-Born-Infeld action, but the resulting 
equations of motion reduce to those derived from (ED).) Since each string is bounded by the two 
junctions, the a parameter runs between two time-dependent bounds 

a q e[s A , q {t),s B , q {t)]. (2) 

The conformal gauge constraints can be written (with ' and • standing for derivatives with 
respect to a and t) as 

X'-Xg = (3) 



q 

i 

and as usual, away from the junctions, the wave-like equation of motion x„ — x" = yields 



< 2 +± q = i (4) 



* q (cr,t) = -(a q ((T + t)+b q (a-t)), (5) 

with &' q = h' q = 1 in order to satisfy the gauge constraints. This resembles the solution for 
standard closed loops. However, one must also take in account the second term of ([T]) which 
imposes, via the Lagrange multipliers f q , that the 3 strings meet at the positions of the junctions 
X^(t) and Xfi(t). Causality requires |Xj(t)| < 1, which together with (TjJ yields [8] 

\SA,g\ < 1 , \s B ,q\ < 1- (6) 

The time evolution of SA, q (t) and SB, q (t) is given by the equations of motion at the junction, 
yielding [8] 

/ii(l - s B ,i) Mi(l - cb.i) (7) 
E q N T, q M q (l-c B>q )> U 

where 

c B ,i = b' 2 (s 2 - t) ■ h'z{s z - t) , Mi = n\ - (jn - /i 3 ) 2 > , (8) 

and similarly by circular permutation. As a result of the presence of the junctions, loops now 
evolve non-periodically in time. 

Like field theory cosmic strings, superstrings can contain cusps — points moving at the 
speed of light |x 9 | = 1 (and hence x' = 0) with x" ^ [lOl |9] - - and kinks, which correspond 
to a discontinuity in x' g . Here we aim to study the gravitational wave emission by such localized 
sources. We first calculate the GW signal in the local wave zone of the source, at distances from 
the source that are larger than the wavelength of interest but smaller than the Hubble radius. 
In this regime we can take the spacetime to be asymptotically flat: g^ u = r]^ v + h^ u , where 
hfj, v «C 1 is the metric perturbation generated by the source. The subsequent propagation of 
the gravity waves on cosmological scales in a Friedmann-Lemaitre spacetime is discussed at the 
end. 
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In a suitable gauge the GW are described by the transverse traceless (TT) part of the linear 
perturbation h ( ^ T) of the spatial metric. This satisfies the linearized Einstein equations, 

Uh% T) = -16ttGT< tt) . (9) 

where T^- is the (TT) part of the stress-energy tensor of the source. In the local wave zone, for 
closed loops of characteristic size L localized around the origin, the solution is given by [6] 

h^fr U ) = ^e^2V T >(um,u>) (10) 

where n = x/||x|| and r = ||x||. Note that h^ T) is actually the time Fourier transform of 
h ( ^ T) and T^J T) the space-time Fourier transform of T^J T) : we will use these abusive notations 
throughout the text. The wave-zone conditions are [TBI [16] 

r > L , r > L 2 oj , r > l/u . 

Hence to obtain the GW emission in this regime it suffices to evaluate T^ T) at points (k, to) 
in the Fourier domain that satisfy the dispersion relation k = urn, where n is a unit vector 
pointing from the source towards the observer. 

Let us put aside the (TT) projection From now on we concentrate on the calculation of 
and merely indicate where the (TT) projection projects out a leading contribution. From the 
action (TjQ) one gets [8] 

r f S B,q(t) 

T^(x ,t )= to dt (x^-xX^I^^Hxo-x.^t))^-^, (11) 

string q s A, q {t) 

namely a sum of contributions of all strings. Each term can be evaluated independently once 
the coupled dynamics of the system has been determined. From now on, we focus on one of 
these contributions and drop the q subscript indicating the string. Moving into Fourier space, 
we have 

r r s B{t) 

T y (k,u;) = n I dt I (x^-x^^e^-^^da (12) 



r r s B(t) 

^ dt a'^a + ^b'^a-ije^-^^da (13) 

2 J JsA(t) 



ss(t) 

For the periodic loops with no junctions considered in [6], changing integration variables to 

u — cr + t , v — a-t. (14) 
and using the periodicity leads to a factorized expression of the form 



Tij oc ( / 4(u)e^«- k - a ("))du)( / b'^e-^+^^dv). (15) 



With junctions, the absence of periodicity prevents one from writing the stress-energy in a 
similar convenient form. Nevertheless, the same change of variables still proves useful to study 
(see figure [5]). Indeed, after changing to (u, v) in (TT3"]) and evaluating at k = am we obtain 



.. poo . ru B {v) 

T ii (w,o;n) = j/ b'^e-^+^^U a^u)e^ n <^du) dv. (16) 



Figure 2: Worldsheet of one of the strings 



Notice that the u and v integrals are no longer factorized due to the v dependence of the 
bounds in the u integral. Indeed, the bounds ua(v) and Ub(v) are defined as follows. Since 
\sa\, \sb\ < 1, for each v there is a unique tA and a unique ts that satisfy sa^a) — tA = v, 
sb^b) ~t B = v. Then, u A {v) = s A (t A ) + t A and u B {v) = s B (t B ) + t B . 
We note for further use that, for J = A or B, 



This means that for the "usual" values of \s'j\ not too close to 1, duj/dv is negative and of order 
1. However, one can also have the limiting values duj/dv = — oo when sj — 1 and duj/dv = 
when sj = —1. Since such values correspond to a junction moving at the speed of light, a 
phenomeneon which will turn out to lead to a gravitational burst, we will consider these more 
carefully below. Note, however, that the equations of motion forbid the strings from shrinking 
at the speed of light (either at junction A: sa = 1 and duA/dv = — oo, or at junction B: Sb = — 1 
and dus/dv = 0) [UJ. Hence one only needs to consider the case of a string expanding at the 
speed of light, again either at junction A with sa = —1 and duA/dv = 0, or at junction B with 
sb — 1 and dus/dv = — oo). 

We are interested in the production of gravitational bursts, that is in the high frequency uj 
regime. We can therefore restrict attention to the u — > oo limit of the integrals in (1161) . where 
standard techniques have been developed to evaluate ffTUj) [HI H3]- We note however that an 
analysis of this kind does not include the low frequency gravitational radiation from the slow 
motion of the string itself, and therefore the resulting stochastic background of gravitational 
waves arising from this. 

3 High frequency behaviour of T XJ : General discussion 

Let us consider integrals of the form 



in the u — ► oo limit. We will use the following standard result (see, for example, [T3]): 



duj sj + 1 



C . 



(17) 



dv Sj — 1 




(18) 
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If 



(i) Vk^ 0, f {k) (a) = / (fe) (6) and (fe) (a) = ^ k) (b) 

(ii) f and are C°° on [a,b] - that is, f(t) and <p(t) are smooth in the interval; and 
(Hi) V£ G [a, b], <f)(t) 7^ - that is, there is no stationary phase or saddle points, 

then I(uS) — > as u — > 00, faster than any power ofl/u. 

(To see this, change the variable of integration to <p using (Hi) and integrate by parts, (i) means 
the boundary term vanishes. Repeating this procedure N times shows that I(oS) = 0(uj~ n ).) 
Thus integrals of the form (]18p are exponentially small for large values of u if the conditions 
(i) — (Hi) hold and the integrand has a rapidly varying phase, i.e. uj ^> \ f/f<j)\. 

Integrals of this type appear in the expression for the stress energy, eq. ( fl5|) and ( TT6l) . In 
the following we will therefore be interested in situations where (at least) one of the above 
conditions does not hold, since these might lead to high frequency contributions to the stress- 
energy proportional to small powers of But we first summarize the results obtained in [B] 
for the factorized stress energy tensor (1151) of cusps and kinks on periodic strings. In this case 
the above results can be applied independently to each integral, and since (ii) is ensured by 
the periodicity, the main contributions to Tjj appear when 

• there is a saddle point in each integral ((Hi) violated in both integrals), leading to a 
contribution oc 1/ u; 4//3 and corresponding to the physical situation of a cusp emitting 
around one specific direction; 

• there is a discontinuity of a\ and a saddle point in the integral over v (or vice versa), 
leading to T^- oc This is the case of a kink emitting in a one-dimensional fan-like 
set of directions throughout its propagation. 

All other cases (discontinuities in higher derivatives of a 1 or b' for example) give smaller contri- 
butions to Tij. 

In the case of loops with junctions we consider here, the stress energy tensor given in eq. ( fT6i) 
cannot be factorized. It can only be written in the form 

T^u^n) = j[°° b' 3 (v)e~^ +n - h ^h(v) dv , (19) 

where ^ 

h(v) ■= [ * <(u)e-^ (u) dn (20) 

Ju A (v) 

with 

#«) = ~(u-n.a(u)) . (21) 

The previous results can only formally be applied to the integral in Ii(v): indeed, its bound 
depend on v, so that the v -integral in (fT9"j) receives contributions that could prevent it from being 
of the same type. However, as we now discuss, in most situations of interest these additional 
contributions have a simple dependence on v, allowing us to generalize the results of [6]. 
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4 High frequency behaviour of T^i first integration 

Consider first the integral h{v) given in (1201) . In order to be in the rapidly varying phase 
regime, we require lu 3> ||a"||. If the loop is not too wiggly, one can assume that ||a"|| ~ 1/L 
(recall that ||a'|| = 1), in which case the condition becomes uj ^> 1/L (note that this is an 
additional constraint on top of the wave- zone conditions (ITTT) ). We now examine the high 
frequency contributions to I%{v) resulting from the violation of one of the conditions listed in 
the previous section^. 

• (i) is violated i.e. sl'(ua(v)) ^ &'(ub(v)). 

This will generically be the case for a string bounded by junctions. Hence this is a novel 
contribution with respect to [6]. After integrating by parts, the leading order contribution 
scales as Co"' 1 and given by 

u=ub(v) 

(22) 

u=ua(v) 

In most cases of interest one has |^- (v) | ~ 1, so that 1 ^a'{uj(v)) * s a s l° wr y varying function of 
v and e~* w ^( Wj ( t ')) is a rapidly varying phase term that must be taken in account when performing 
the integral over v in (]19j) . though the latter remains of the same type as (|18[) [E 

• (ii) is violated in the following way: there exists a e]ua(v), ub(v)[ where a[ is discontin- 
uous. (We could also consider non regularities in higher order derivatives of a, but they would 
lead to higher order powers in 1/uj.) 

To evaluate this term, we first rewrite 

k{v) = I * a! i {u)e ibJ ' t> du + / * a'^e^du, 

J ua(v) J u* 

and then integrate once by parts as above to find 

I disc {v) = — | a 'i( U *} a i( U * ) | e iuj(f>(u,) ^3) 

ko \ 1 — n ■ a!(ut) 1 — n ■ a!(u7 ) J 

As above, we find an u~ l fall-off. The only dependence on v comes from the bound ua{v) < 
u* < Ub{v). This means that when we will perform the integral over v in ( 1191) in section V 
below, we will need to restrict the v domain of integration appropriately. 

• (Hi) is violated: there exists^ a u s E [ua{v),Ub{v)) where the phase (p(u) has a vanishing 
derivative <fi'{u s ) = —(1 — n ■ a!(u s ))/2 = 0, or in other words n = a!(u s ). 

An analysis of situations where several conditions are violated is beyond the scope of this work. 
2 As mentioned above, when the string expands at the speed of light at junction B, diverges. One might 
expect this leads to a discontinuity-like contribution in the second integral. However, both the derivative with 
respect to v of the amplitude i_°^a'*(u£(v)) an< ^ °^ ^ ae ph ase i n the exponential diverge. The condition tpto ^> ^ 
still holds since cancels on both sides. Therefore, no discontinuity-like contribution is introduced in the 
second integral 

3 Strictly speaking, u s could lie outside [ua(v),ub(v)], but close to ua or ub and still give rise to a non 
negligible contribution (cf Appendix A). 



rboundary . 



too 



1 — n ■ a!(u) 



-iu)(f>(u) 
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In this case the leading contribution to h(v) comes from the vicinity of the saddle point and 
is obtained by Taylor expanding <p(u) and a'^u) around u s . 

To do so, notice that the gauge conditions (on the worldsheet) enforce 4>"(u s ) = 0. Indeed, 
4>"(u) = • a"(u) and since a' (it) 2 = 1, it follows that a' • a"(u) = 0. Furthermore, 4>'"(u s ) ^ 



,"l 1 2 



so that 



because upon taking the derivative of a' • a"(w) = 0, one gets a' • a"'(u) = 
4>'"(u) = — ||a"|| 2 /2. Regarding the Taylor expansion of a'^u), the first term is a'^Ug) = rii 
which is suppressed by the (TT) projection operator (it amounts to a gauge term [6]). To 
summarize, we therefore need to Taylor-expand the phase to the third order and a'^u) to the 
first order: 

(24) 

a'^u) ~ a'^Us) + {u- u s )a'l(u s ) , 
Thus, up to a gauge term proportional to rij, 



u B (v) 
u A (v) 



h ~ e -^ w / (u - u s )a'^u s )e-^"' M{n) du 



~ -e-^ Us) a'!{u s 



we~ m ' dw 



2/3 



a'l(u s )e- iuJ ^ Us) 




we~ lw dw 



w 



u\r(u s )\\ 1/3 



\u - u. 



(25) 



where we have used f°° we %w dw = 
third line, we have assumed that the saddle point lies far from the boundaries: 



1 T(|). We note that in going from the second to the 



\u s -uj[v)\\^"[u s )\^->\ 



(26) 



thus allowing the domain of integration to be extended from — oo to +oo. This will not always 
be the case, however, and the general result reads 



r saddle , 



v) = gauge term + 



^-<K)e^- n ' aK)) 



12 \ 2/3 / i _,2 



n ■ a" 



^rig) ) C(W A ,W B ) 



with 



(27) 



C(w A ,w B ) 



w e %w dw 



J: 



w e~ 



" 3 dw 



B{w A )-B(w B ), B{wj) 



rwj 
J — oo 



w e tw dw 



we™ dw 



(28) 



and 



w A (v) 



uj\n ■ a 
12~ 



1/3 



(u s -u A (v)), w B (v) 



uj\n ■ a 



1/3 



(u B (v) - u & 



The behavior of the function C(w A ,w B ) is studied in detail in Appendix A. Its important 
features are the following, 
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• when the saddle point lies well in between the upper and lower bounds of the u-integral 
C(wa(v),wb(v)) obviously reduces to 1, as in (1251) . 

• when the saddle point is located far outside the domain of integration then C(wa(v),wb(v)) 
reduces to zero. 

• when the saddle point is near ua(v) (resp Ub(v)), wa(v) (resp wb(v)) is of order one. 
C(wa(v),Wb(v)) then reduces to B(wa{v)) (resp B(wb(v))) given in figure El We note 
that B{wa) — > 1 for large values of wa, with the envelope of its oscillations decreasing as 
1/w. 




Figure 3: Real and Imaginary parts of B(wa) as a function of wa 

The entire dependence on v in ( 1271) is contained in C(v) = C(wa(v),wb(v)), which is a 
smoothed version of the step function 6{u s — ua(v))9(ub(v) — u s ) (cf Appendix A). This enters 
the integral over v as follows. Let va, s (resp. vb, s ) be the value of v for which ua(v) = u s (resp. 
ub(v) = u s ). Since ua and ub are decreasing functions of v, C(wa(v),Wb(v)) is actually a 
smoothed version of the step function 6{v — va, s )@(vb, s — v). But we ought to verify whether 
C is rapidly varying (in the vicinity of Va, s and Vb, s ) compared to the phase appearing in 
e -f («+n.b(«))_ In ty p ica i cases w h ere | « l,|^f | « 1, one has « ~ O(l) in 



the region of interest (see Figure [3]). Hence 

^> « ^, s ) « ^(«, « ^)) (^|0'"(n s )|) 1/3 ^ « (^|0'"(n s )|) 1/3 (29) 



For strings that are not too wiggly we expect | <p'" | ~ | |a" 1 1 ~ 1/Zr so that 

(v « u A>a ) « — < (30) 



since we assumed uj 3> Hence C(v) enters as a multiplicative, slowly varying amplitude 

in the integral over i@. In particular it does not introduce a boundary term in the integral 
over v (unless of course dus/dv diverges, see below), which implies junctions do not radiate 
spontaneously. 



4 Strictly speaking, we need to prove that («) is slowly varying in the sense that (I (v)b'j (v)) jj^rj^ *C 

uj. After applying the derivative one is left with two contributions: the first one is ^ which, as before for the 
a' term, is of order <C u>. The second one is -j = ^ which, according to what we just said, is -C to. 
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5 Full high frequency behaviour of XL 



We are now in a position to perform the remaining integral in (fT9l to obtain the high frequency 
behavior of the various contributions to the stress energy source of GW bursts. The analysis of 
the previous section shows that, in all cases, the integral over v is of exactly the same type as 
the integral Ii(v) over u. We can therefore apply the same methods. Generalizing the definition 
we used above in the saddle point case we denote by va{u) the value of v for which ua(v) = u 
(and the dual definition for B). Since we are interested in the high frequency regime, we list 
the different contributions with increasing powers of 1/uj. 

5.1 Contributions in l/o; 4 / 3 

• saddle point in the u integral at u s / saddle point in the v integral at v s (standard cusp) 

The contribution from the integral over u can be treated as a slowly varying amplitude. 
Hence the integral over v is formally of the same type as the it-integral, with the product 
C(v)bj(v) now acting as the slowly varying amplitude. The derivative of this evaluated 
at v s contains a term proportional to b'j(v s ), which is a gauge term for the gravitational 
waves, and a term proportional to C(v 3 )bfj(v s ). Our final result reads 



Tij « -^ a J / ( Us )6j(i; a ) e T[««-''.-»-W".)+b(«.))] 

/ 12 \2/3/ 12 \2/3 1 O, , 

(31) 



n = a!(u v ) = —b'(v 



a!"(u s )\\ 1/3 , fu\n-a!"(u s )^ 1/3 



with 
and 

wa = ( — i2^ Sn ) ^ Us ~ UA ( Va ))> w b = ~ ( — i2^ Sn ) ( u b(v s ) - u s ). 

This is the case of a cusp which emits in the direction given by a!(u s ). If the cusp occurs 
away from the junction, C reduces to 1 and we recover the standard result of Damour 
and Vilenkin [6]. 

5.2 Contributions in 1/cj 5//3 

This is the first case in which we find novel contributions specific to strings with junctions. 
• discontinuity in d i at some / saddle point in the v integral at v s (standard kink) 



This is the standard case of a left-moving kink propagating on the string and emitting in 
the fan of directions n = —h'(v) generated by the right moving waves: 




12 \ 2 / 3 



2v 3 3 
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(32) 



with 

n = -h'(v s ) 

and 

u>a=( ^ J ( u » - u a (««))> % = - J (^bK)-^). 

In (j3~2~|) . C is the complex conjugate of C (induced from the fact that in the Taylor 
expansion of the phase of the v integral, we now have a positive 4>'"(v), and hence this is 
equivalent to changing % to —i in all integrals.) 

In order to illustrate how the presence of the function C(wa,wb) modifies the case of a 
standard loop with no junction, let us consider a left moving kink propagating between B 
and A and ask in which directions n (described by points on the 3D unit sphere) it emits. 
The set of points {—h'(v)} draws a curve C on the Kibble- Turok sphere, as illustrated in 
figure HI Away from C, Ty(n) ~ whereas on the curve the amplitude is given by (13*2*1) . 
The interval between positions va(u*) and Vb{u*) on the curve contains the right moving 
waves —h'(v) effectively seen by the kink while it propagates, so we expect to have an 
emission only in or around these directions. On this interval - and not too close to the 
endpoints - we have C ~ 1. Far from this interval (in blue on fig HJ, C ~ 0. Finally, 
C smoothly decreases from 1 to in a short interval Av of a few (L 2 /^;) 1 / 3 near va(u*) 
and vb{u*). Hence the overall probability that a GW burst comes with a value of C 
significantly different from one is of the order (La;) -1 / 3 <C 1. We conclude therefore that 
up to small corrections the predictions of [6] apply to GW bursts from kinks on cosmic 
strings with junctions. 




Figure 4: Directions of emission of a propagating kink on the Kibble- Turok sphere 



• saddle point in the u integral at u s j discontinuity of b^- at some (standard kink) 



5 We note that although this is not taken into account by (f5!""|) . smoothly vanishes when n moves away 
from C. 
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This is a situation that mirrors the previous one, namely a right-moving kink propagating 
on the string and emitting in a fan-like set of directions given by n = a'. The formulae 
are similar to those above. 

• boundary term (A) in the integral over u / saddle point in the second integral at v s 



Inserting ( 1221) into ( 1191) . we obtain: 



rr, f x ¥ f°° W a 'i(M*>)) -t 
1,-Au, uii) — — / — - ; — e 



(v+n-b(v)) - (ua(v)— n-a(uji(«)) 



dv (33) 



wnere i-n-a'(MA(^)) * s a S ^ 0W ^ varying amplitude compared to the rapidly varying phase 



<t>( v ) = 2 i v + n ' b ( v )) ~ ( u a( v ) ~ n • sl(ua(v)) 



The saddle point condition is that 

1 



duA 



1 + n • b'(v a ) + -t^K) (n • a'(us(v a )) - 1) 



0. 



(34) 



(35) 



Since both terms in the sum are positive, they both need to vanish at v s . Moreover 
n ■ a!(uA{v s )) 7^ ~ 1, otherwise we would have a cusp at the junction. Hence we are led to 
the following conditions: 



n • b'(v 



n 



-b'(v x 



dv 







(36) 



The second condition means the string is expanding at the speed of light at junction 
A. The other case where the string expands at the speed of light at junction B would 
correspond to ^p- = — oo, and thus to a discontinuity in the v integral. We will treat this 
separately below. 

The amplitude of this contribution can be calculated using the saddle point treatment 
described previously (namely a Taylor expansion of the amplitude and the phase). We 
have 




^(v s ) (n • a!(u A (v.)) - 1) + n • a>») (^(v.) 



=o 



— 

=0 



(37) 

The first term vanishes because of the gauge conditions on the worldsheet, and the last 
term vanishes due to the second condition in (|3"B"]) . The second term also vanishes. Indeed, 
since -^-(v) ^ and ^(1;,) = 0, ^f(f) reaches its maximum at v s . Hence we need to 
Taylor expand the phase to third order: 



n ■ b"'(v s ) + («.)(* ■ a!(u A (v.)) - l) 



(38) 
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On the other hand, Taylor expansion of the amplitude to first order gives 

a'^UAjv)) +( v - v )h"(v ) giM^lj 

W 1 - n ■ a'(u A (v)) gaUge t6rm + (U Vs)bj {Vs) 1 - n ■ a'(tu(«.)) 

where the gauge term is proportional to b'j(v s ) = —rij. 

Hence the leading term of the energy-momentum tensor is given by 

Tij(u), urn) « —b^vs)- — f-^rW {v s ) ) 

uj b i A J 1 — n ■ &'(ua{v s )) V / 



(39) 



6 2/3 o 

2v/3 V 



(u s +n-b(« s )J - (uA(t's)-n-a(uA(fs)) 



(40) 



This contribution corresponds to the situation in which the string expands at the speed 
of light at junction A and emits in a direction corresponding to n = — h'(v s ), with an 
amplitude proportional to 1/u; 5 / 3 . We now turn to the case where the string expands at 
the speed of light at junction B. 

saddle point in the u integral at u s / string expanding at the speed of light at junction B 



In this case dus/dv diverges, which gives rise to a discontinuity-like contribution in the 
integral over v. From ffT§]) and (1271) we have 

j 4^/3 a *^Je {\n.a!"(u s )\) V 

C(v)b'Av)e-^ iv+n - h{v)) dv , 



where C(v) = C(wa(v), wb(v)) was defined in f[2"H]) . We are interested in cases in which a 
discontinuity in C (or at least a variation faster than that of the phase in the exponential, 
\dC/dv\ 3> oj) is induced by the divergence of dus/dv or equivalently duus/dv at some v*. 
Yet, to achieve such a discontinuity in C, we need v* to satisfy 

KWI<1. (42) 

Indeed, in this region we have OC/Owb ~ 1 whereas elsewhere C does not vary with 
wb at leading order so that a strong variation of wb{v) would have no effect on C(v). 
Physically, this requirement means that the saddle point condition is satisfied near the 
junction at a moment when sb reaches 1. 

From now on we assume that (1421) is satisfied. For \dC/dv\ 3> u, we can treat C as being 
discontinuous at v* where dC/dv diverges. In this case, the leading contribution is the 
same as that of a discontinuous C with a gap AC(v*) equal to the variation of the true 
function C over the interval where \dC/dv\ ^> oj: 



(43) 



2^3 v 3 
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We therefore need to evaluate the order of magnitude of AC(v*). Let us start by de- 
termining the interval around over which \dC/dv\ ^> uj. As was already explained in 
(12TI1) and flHUl) . \dC/dv\ w Irfn^/d^Kcj/L 2 ) 1 / 3 so the condition |dC/rfv| 3> translates into 
\duB/dv\ 3> (a>L) 2 / 3 . Since I^mb/g^I = (1 + s B )/(l — s B ) w 2/(1 - sb) (because ~ 1 
in the region of interest), this finally yields the condition 

l - iB<< !^- (44) 

We therefore need to determine how sb(£) behaves around a point t* where suit*) = 1- 
To this end we investigate the dynamics of a junction using the equations recalled in 
Section [5] (we momentarily drop the B index in equation ((7|), and assume the above 
discussion applies to string 1, namely si(t*) = 1, while strings 2 and 3 are assumed not 
to be expanding at the speed of light.) For convenience define t>2* = s 2 (U) — t* and 
'^3* = s 3 (t*) — t*. Then, from ([7]), ci(£*) = 1 so that b 2 (t>2*) = b 3 (u 3 *) = B. Expansion of 
b 2 (t> 2 ) and b 3 (t>3) around v 2 * and t> 3 * respectively then yields 

b' 2 (v 2 ) « B + b' 2 '(^)(t; 2 -^)+b' 2 "(^)^ 2 



b> 3 ) ~ B + b 3 > 3 *)(?; 3 -^+b 3 > 3 *) 



2 

(«3 - ^3*) 2 



where we have used the fact that b^^*) • B = h'^(v^) • B = because of the gauge 
constraints on the worldsheets. Thus 



ci{t) ~ 1+ 2 b 3 ^ 3 *) ' B + 2 2 ^ ^ ' (45) 

+ (s 2 (t) - t- v 2 *)(s 3 (t) - t- v 3 *)b 2 (v 2 *) ■ b 3 (?; 3 *) 

We can also expand s 2 (t)—t ps U2*+(s2(**) _ !)(*— **) an d s 3 (t)-t ps u 3 *+(s 3 (^)-l)(t-^) 
where generically (s 2 (t*) — 1) and (s 3 (t*) — 1) are of order 1. Then, using the fact that 
the three scalar products involved in (1451) are of order 1/L 2 , we can write, around t* 

ci(t)-l^--^(t-t*) 2 (46) 

where the minus sign ensures that c±(t) ^ 1. Plugging this into equation (JTj), we finally 
get the expansion around t* 

l-aiW*-!^-*.) 2 - (47) 

Therefore, (l44l) holds over an interval of time around t* of size At of order At pa 
L/^luL) 1 / 3 = L 2 / 3 c<j _1 / 3 . Since «s = s#(t) + t, the variation of ub during this inter- 
val At is of the same order as At i.e. Aub ~ L 2 / 3 ^;" 1 / 3 and so the variation of wb is 
Aw B ~ (u;L~ 2 ) 1//3 Am b p^ (cjL" 2 ) 1 / 3 /, 2 / 3 ^" 1 / 3 ^ 1. Finally, since we assumed that w B (v*) 
is in the region where OC/Owb ~ 1, this yields 

ACO*) ^ 1 (48) 



To conclude, if n = a!(u s ) and dus/dv diverges at some v* such that 
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\w B {v* 



( ^ n a i2^' s ^ ) \ u b( v *)~ u s\ ^ 1 then the integral over v receives a discontinuity- 



like contribution leading to 

[A ... . b'Jv*) ( 12 \2/3 

(49) 



' 1 + n • b' 



x 



— Lr(_) x Ac(^) e ^( Ws -^- n - (a(Us)+b ^ )) ) 

2v 3 3 



where with respect to (1431) we now write AC(t>*) on the second line to indicate that this 
is a numerical factor of order 1. The string expands at the speed of light at junction B. 



5.3 Contributions in 1/u 2 

• Discontinuity in a' { at some it* / Discontinuity in b'^ at some v* G [va_{u*),vb{u*)] 



T 



%0 u 2 \ 1 — n • a!(ut) 1 — n ■ a!(u* 

(50) 

ty( V * ) ^'A V * ) \ e ^[«*-«*-n-(a(«,)+b(t;,))] 

1 + n-b>+) 1 + n-b'(y-) y 

This is the situation where a left moving and a right moving kink meet during their prop- 
agation on the string and emit in all directions. 

Discontinuity in a' { at some / boundary term in the integral over v 



rp _ _£_ ( <«) <K ) 



uj 2 \ 1 — n ■ a'(u^) 1 — n ■ a'(it* ) 

/ x 7 (51) 

/ b'^VB^)) \ ^ lUt _ VBM _ n . (a(Ut)+h{VB(Ut))] 

^1 + ii'b'Ma,)) J 

This corresponds to a left-moving kink passing through a junction and emitting in all 
directions in space. 

boundary term (A or B) in the integral over u /discontinuity of b'j at some 



x I - - -- - ■ ■ 1 glf (v*)—v* —n-(a(uA (v*))+h(v* )] 



uj 2 \l + n • b'(u+) 1 + n • b'K 

(52) 

a-(M^)) 
1 - n • a'(uA(v*)) 

This describes a right moving kink passing through a junction (A here) and emitting in 
all directions in space. As one expects, the symmetry with respect to the case of a left 
moving kink passing through B is recovered despite the explicit breaking of symmetry 
made by choosing to integrate on u first. 
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6 Conclusion 



The transverse traceless projections of the various contributions to the stress-energy tensor 
found in Section V fully determine the GW emission from cusps and kinks in the local wave 
zone of the source through eq ( ITUi) . The observed signal is obtained by parallel propagation 
of the gravity waves in a cosmological background, along the null geodesic followed by the 
GW. This gives rise to the usual redshifting of time intervals between emission and reception, 
which in the Fourier domain corresponds to f em = (1 + z)f rec . One obtains the following 
order-of-magnitude estimates for the logarithmic Fourier transform of the observed amplitude 
of individual GW bursts emanating from the various high frequency sources discussed here 

W~> ~ ((1+^5= ~^ m 

where to denotes the present age of the universe. Most importantly, the exponent a in this 
expression is determined by the high frequency behavior ~ co>~( 1+a ) of the stress-energy sources 
calculated above. 

In the time domain fl53|) corresponds to a signal of the form ~ \t — t c \ a , where t c is the 
arrival time of the center of the burst. The total GW signal of a superstring loop consists of the 
sum of the various bursts and a slowly varying component due to the low frequency modes of 
the strings. Despite its vanishing at t = t c , bursts are distinguishable from the slowly varying 
component because the curvature associated with ( l53i) diverges as ~ \t — t c \ a ~ 2 , exhibiting 
clearly the spiky nature of GW burst^l 

We are now in a position to summarize the modifications of the GW signal from bursts 
induced by the presence of junctions, with respect to the case of standard loops. As expected, 
away from the junctions we recover the same GW signal as Damour and Vilenkin (6j [7j, both 
for a kink propagating on a string - for which a = 2/3 and the emission occurs in a fan-like 
set of directions - as well as for a cusp, where a = 1/3 and the emission occurs in a small cone 
around a particular direction. 

However for bursts from sources near junctions we find the predicted GW amplitude is 
modified by a smooth correction factor C. This means that for GW bursts from cusps and 
kinks near junctions one effectively measures a combination of the tension \x and C. Since these 
are rather improbable events, however, the overall effect of this on the (statistical) predictions 
of bursts from cosmic superstrings is likely to be small. The smoothness of the correction factor 
also implies that despite their spiky nature, junctions do not radiate spontaneously (at least 
for strings without structure, as considered here). 

Furthermore, the presence of junctions gives rise to several novel sources of GW bursts, 
most notably from 

• a string expanding at the speed of light at the level of a junction, which emits in a specific 
direction with an amplitude (1531) with a = 2/3. 

• a kink passing through a junction, which emits in all directions with an amplitude (1531) 
with a — 1. 

Finally, the situation where a left-moving and a right-moving kink pass through each other 
during their propagation remains unchanged and has a — 1. 



6 In practice, the observer will never lie exactly in the direction of emission. However the signal is detectable 
inside a small cone around this direction. This introduces a high frequency cut-off on the waveform or, equiva- 
lently, a smoothing of the signal around t c in the time domain. 
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To translate our results into observable waveforms one ought to sum the individual contri- 
butions from cusps and strings in a cosmological network of string loops. However the rate of 
occurrence of events of different types strongly depends on the dynamics of the loops and on 
the evolution of the network, which is currently poorly understood. This gives rise to significant 
uncertainties in the resulting GW signal. Indeed, an increase in the number of events of a given 
type lowers the threshold redshift for observation and hence reduces the expected dilution of the 
signal. The individual contribution with the smallest value of a therefore need not necessarily 
provide the dominant contribution to the observed signal. 

In particular, concerning the case at hand, one might expect junctions to enhance the 
number density of kinks whereas it appears unlikely that the rate of occurrence of cusps is 
significantly affected by the presence of junctions [12]. This would mean that even though 
individual GW bursts from cusps are stronger, the increased number density of kinks on loops 
with junctions might compensate for the difference in strength. Since we have identified GW 
bursts from events involving kinks that are specific to strings with junctions, this would provide 
a promising route to observationally distinguish between gauge theory cosmic strings with no 
junctions and superstrings. 
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A Detailed study of the function C(v) = C(wa,wb 

In (1271) . we gave the general result for the integral ( 1201) due to a saddle point: 



T saddle i 



gauge term + 



2/3 



12 



n • a" 



2/3 



C(w A ,w B ) 



where 



C(w A ,w B ) 



WA w e~ iw dw 

WB 



I: 



w e 



dw 



B(w A )-B(w B ), B(wj) 



WJ w e~ iw ' 

oo 



dw 



w e~ 



dw 



and 



w A iv) 



u\n ■ a 
12~ 



1/3 



(u s - u A (vj), w B (v) = - 



uj n • a 
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1/3 



(u B (v) - U s ). 



(54) 



(55) 



(56) 



The function C depends on v because the position of the saddle point relative to the bounds 
u A (v) and u B {v) of the integral does. This Appendix contains a detailed analysis of C(v), 
whose behavior needs to be understood in order to determine the saddle point contributions in 
the three cases listed under f j28|) . 



Envelope of oscillations of B 

When the saddle point is far from w B , we can write C(w A , w B ) ~ B(w A ), where B{w A ) is 
displayed in Figure [3j We first show that the envelope of the oscillations of B scales as l/w A . 
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Consider for example the real part $1(B(wa)) = ol f_ A w sm(w 3 )dw, where 
a = — i/{ we~ %w3 dvo), which oscillates around its asymptotic value of 1 with decreasing 
amplitude. Let Wk = (nk) 1 ^ 3 be the values where its derivative (aw sm(w 3 )) vanishes. These 
points correspond to the relative maxima (for odd values of k) and minima (for even values of 
k) of 31(B(wa))- Then, because Wk+i — Wk is small, the amplitude of the oscillations around a 
certain Wk (say k even for example) is given (for large values of k) by 



B(w h+1 ) - B(w k ) a ■ , 3m « 



2 / 2 



cos [w 



wsin(w )dw ~ — ~ - — . (57) 



3w 



3w k 



We infer that B(w) is of order 1/w in the limit w — > — oo. In particular, if wb is negative 
and large in absolute value, we have C(wa,wb) ~ B(wa)- 

Contributions to ( 127)) due to a saddle point inside [ua(v),ub(v)} 
Following the discussion in section III, we investigate the contribution due to a saddle point 
lying between ua{v) and ub(v). This corresponds to studying C(v) between the points Va, s 
and vb, s defined at the end of Section IV. Since 4>'"(u s ) ^ 0, u s e [ua{v),ub(v)] translates into 
wa{v) ^ and wb(v) ^ 0. We need to distinguish between two regimes: 

a. Saddle point far from both bounds: \u s — mj(v)||o;0" / (m s )| 1 / 3 3> 1 

Then we have wa{v) 3> 1 and (v) | ^> 1 so according to the discussion above, C(v) ~ 1 
up to corrections of the order of 1/w a, 1/ujb <C 1 

b. Saddle point near one of the boundaries (e.g. A): \u s — ua{v)\\u}())'" '(us)] 1 ^ 3 < 1 

First, note that the saddle point cannot be close to both bounds at the same time: (uxj)'") 1 ^ 3 
is of the order of uL~ 2 3> l/L so \u s — Ua(v)\ L; since \ua(v) — Ub(v)\ is of the order of L, 
\u s — Ub(v)\ has to be of the order of L too so we still have |wb(?;)| ^> 1. 

Therefore, up to a negligible correction of the order of 1/wb, C(wa{v), wb(v)) ~ B(wa{v)). 
This time wa is of order 1 (and still positive). The behaviour of B(wa) is plotted in Figure El 
The contribution in this case is of the same order as the one far from the bounds but the result 
is multiplied by the complex factor B(wa) of order 1. 

Contributions to ( 1271) due to a saddle point outside [ua(v),ub(v)] 
Up to now we calculated the leading contribution from a saddle point present anywhere 
inside the interval [ua{v),ub{v)]. Does this mean that for the values of v for which the saddle 
point lies outside this interval, I{v) discontinuously drops to zero? 

We argued in Section [H] that, if Uiv) obeys all conditions (i)-(iii), then it is negligible 
because it vanishes exponentially in the infinite u limit. However, since u is actually large 
but finite, there might also exist non negligible contributions if I(v) is only close to satisfying 
these conditions. In particular, if some saddle point lies just outside the interval of integration, 
the derivative at the nearest bound is very close to vanishing so we can expect to have some 
contribution. Another way to say this is that the contribution comes from an interval around 
the saddle point of size ~ (L 2 /u) 1 ' 3 which can intersect [ua(v),ub(v)] even if the saddle point 
is outside it. Obviously, in the case where the saddle point is very far from the interval of 
integration, the contribution needs to be negligible. 

a. Saddle point outside [ua(v),Ub(v)] and near one of the boundaries (e.g. A): 

\u s — UA{v)\\uj(t)"' (u s )\ l l 3 < 1. Here, wb ^ 0. For the same reasons as before, we have \wb\ ^> 1 
and C(wa(v),wb(v)) ~ B(wa(v)). The only difference is that now wa ^ 0. This is why in 
Figure [3] we also plotted B for negative values of its argument. 

b. Saddle point far from [ua(v),ub(v)] 
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In this case, both wa and wb have the same sign and their absolute values are ^> 1. Accord- 
ing to the discussion in the first section of this appendix, we have at leading order C{v) ~ up 
to corrections in 1/wa, 1/Vb <C 1. 

Summary 

It should now be clear that C(v) is a smoothed version of 6{v — va, s )9(vb, s — v). In Figure 
[5] we sum up the different cases and show the typical behaviour of the function C{y) provided 
that the derivatives ^ are of order 1 which is the generic case. 




Figure 5: Typical behaviour of the function C(v). The different cases are represented on the 
left panel. The left and right brackets respectively indicate the position of Ua(v) and ub(v) on 
the u axis while the brown tick is the position of u s 
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